Primary exercises
- Dietary intakes. (Create a vector, use it in
calculation.)
Four patients had daily dietary intakes of 2314, 2178, 1922, 2004
kcal.
Make a vector intakesKCal of these four values.
What is the class of this vector?
Convert the values into in kJ using 1 kcal = 4.184 kJ.
intakesKCal <- c( 2314, 2178, 1922, 2004 )
intakesKCal
[1] 2314 2178 1922 2004
class( intakesKCal )
[1] "numeric"
intakesKCal * 4.184
[1] 9681.776 9112.752 8041.648 8384.736
- More dietary intakes. (Combining/appending/merging
vectors.)
Additional set of intakes is provided: 2122, 2616, NA, 1771 kcal.
Use c() to append the new intakes after values in
intakesKCal and store the result in
allIntakesKCal.
Print the combined vector and print its calculated
length.
intakesKCal2 <- c( 2122, 2616, NA, 1771 )
allIntakesKCal <- c( intakesKCal, intakesKCal2 )
allIntakesKCal
[1] 2314 2178 1922 2004 2122 2616 NA 1771
length( allIntakesKCal )
[1] 8
- The average and total intakes. (Calculating means and sums,
skipping missing values.)
Calculate mean intake for patients in vector
intakesKCal.
Next, calculate mean intake for patients in vector
allIntakesKCal.
Can you explain the result?
Check help for ?mean, in particular the na.rm
argument.
Use the extra argument na.rm=TRUE to calculate the
mean of non-NA elements of
allIntakesKCal.
Check help for ?sum how to omit NA elements in
sum calculation.
Now, calculate the total sum of allIntakesKCal
intakes ignoring the NA element.
mean( intakesKCal )
[1] 2104.5
mean( allIntakesKCal )
[1] NA
# since one element is missing, the mean is unknown
# ?mean, adding argument na.rm=TRUE will omit NA elements
mean( allIntakesKCal, na.rm = TRUE )
[1] 2132.429
# ?sum also allows na.rm=TRUE argument to skip NA elements
sum( allIntakesKCal, na.rm = TRUE )
[1] 14927
- Selecting valid intakes. (Selecting non-missing elements;
logical vectors.)
Understand the result of is.na( allIntakesKCal ).
Now, negate the above result with ! operator.
Use above vectors as argument to sum to calculate the
number of missing and non-missing elements in
allIntakesKCal.
Understand allIntakesKCal[ !is.na( allIntakesKCal ) ].
is.na( allIntakesKCal ) # TRUE marks positions with missing data
[1] FALSE FALSE FALSE FALSE FALSE FALSE TRUE FALSE
!is.na( allIntakesKCal ) # TRUE marks positions with available data
[1] TRUE TRUE TRUE TRUE TRUE TRUE FALSE TRUE
sum( is.na( allIntakesKCal ) ) # number of missing elements
[1] 1
sum( !is.na( allIntakesKCal ) ) # number of non-missing elements
[1] 7
allIntakesKCal[ !is.na( allIntakesKCal ) ] # keeps elements which are not NA
[1] 2314 2178 1922 2004 2122 2616 1771
sum( allIntakesKCal[ !is.na( allIntakesKCal ) ] ) # same as sum( allIntakesKCal, na.rm = TRUE )
[1] 14927
- Generating random kcal intakes. (Generating normally distributed
random numbers; descriptive statistics.)
The code v <- rnorm( 10 ) would sample 10 numbers from
the normal distribution and store them as a vector in
v.
Print v. Then repeat v <- rnorm( 10 ) and
print v again. Has v changed?
Next, read the manual of rnorm and find how to generate
random numbers with given mean and standard deviation
(sd).
Now, in v simulate kcal intake by generating 15 random
numbers with mean=2000 and sd=300.
Print v and find by eye the smallest and the largest of
these numbers.
Try to use the functions min and max on
v – have you found the same numbers by eye?
Calculate the mean, median and the standard
deviation (sd) of v.
v <- rnorm( 10 ) # a vector of random numbers
v
[1] -0.82446256 1.10541227 0.37879338 -0.04248810 -0.90365982 0.30283131 1.99623212 0.04544249 0.45141852 -0.92189770
v <- rnorm( 10 ) # another vector of random numbers
v
[1] -0.7170064132 -0.5644114018 -1.1892413312 -0.6296518231 1.5734457672 -1.9331310668 -1.1680287880 0.0009461918 0.7974261152 -0.1574395731
v <- rnorm( n = 15, mean = 2000, sd = 300 )
v
[1] 2569.448 2070.651 2095.254 1844.532 2106.756 1745.908 1508.779 2548.740 1289.149 2282.688 1687.922 2331.825 2063.216 2282.436 1818.796
min( v )
[1] 1289.149
max( v )
[1] 2569.448
mean( v ) # is it close to 2000? try several random v vectors and see the effect of growing n
[1] 2016.407
median( v )
[1] 2070.651
sd( v ) # is it close to 300? try several random v vectors and see the effect of growing n
[1] 367.0505
- Selecting and counting some kcal intakes. (Selecting elements by
a condition; logical vectors.)
In v simulate kcal intake by generating 15 random numbers
with mean=2000 and sd=300.
Type v < 2000 and understand the result.
How to interpret the number produced by
sum( v < 2000 )?
How to interpret the number produced by
sum( !( v < 2000 ) )?
v <- rnorm( n = 15, mean = 2000, sd = 300 )
v
[1] 1923.3155 2042.5430 2013.6030 1452.1621 2269.4636 999.6216 1962.1153 2020.9400 2093.6657 2083.8385 1712.6447 1796.3014 2150.5820 1612.0690 1486.3805
v < 2000 # TRUE corresponds to elements of vector v SMALLER THAN 2000
[1] TRUE FALSE FALSE TRUE FALSE TRUE TRUE FALSE FALSE FALSE TRUE TRUE FALSE TRUE TRUE
v[ v < 2000 ] # selected elements of v smaller than 2000
[1] 1923.3155 1452.1621 999.6216 1962.1153 1712.6447 1796.3014 1612.0690 1486.3805
sum( v < 2000 ) # number of elements in vector v smaller than 2000
[1] 8
sum( !( v < 2000 ) ) # number of elements in vector v GREATER OR EQUAL than 2000
[1] 7
sum( v >= 2000 ) # same as above
[1] 7
- Head and tail.
Often there is a need to quickly look at the beginning
(head) or at the end (tail) of a vector.
Try these functions to show the first 5 and the last 7 elements of a
randomly generated vector v <- rnorm( 20 ).
v <- rnorm( 20 )
v
[1] 0.45253519 0.63380039 -1.21741887 0.31068980 0.46470936 -2.11669339 -1.49189600 1.91520679 0.08976667 -0.63424471 -0.86648686 -0.32525920 0.93328004
[14] 0.75450744 -0.98147268 -0.57422432 -0.33942420 -0.08137151 -0.30676734 -0.18507061
head( v, 5 )
[1] 0.4525352 0.6338004 -1.2174189 0.3106898 0.4647094
tail( v, 7 )
[1] 0.75450744 -0.98147268 -0.57422432 -0.33942420 -0.08137151 -0.30676734 -0.18507061
- Elements of a vector.
Let’s assume that eight persons had caloric intakes of 2122, 2616, NA,
1771, 2314, 2178, 1922, 2004 kcal.
Make a vector intakesKCal of these eight values (in the
given order).
Use the square brackets to get the 4th element of
intakesKCal.
Use the square brackets and the colon operator (:) to get
the elements from the second to the fifth (inclusive).
Define another vector poses with values 1, 3, 5, 7. Use it
get the 1st, 3rd, 5th and 7th element of intakesKCal.
Finally, get the 1st, 3rd, 5th and 7th element of
intakesKCal typing numbers directly inside
[...] (without using an extra poses
variable).
intakesKCal <- c( 2122, 2616, NA, 1771, 2314, 2178, 1922, 2004 )
intakesKCal
[1] 2122 2616 NA 1771 2314 2178 1922 2004
intakesKCal[ 4 ]
[1] 1771
intakesKCal[ 2:5 ]
[1] 2616 NA 1771 2314
poses <- c(1,3,5,7)
intakesKCal[ poses ]
[1] 2122 NA 2314 1922
intakesKCal[ c(1,3,5,7) ]
[1] 2122 NA 2314 1922